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The spread of infectious diseases near the epidemic threshold is investigated. Scaling laws for the 
size and the duration of outbreaks originating from a single infected individual in a large susceptible 
population are obtained. The maximal size of an outbreak n* scales as N'^^^ with N the population 
size. This scaling law implies that the average outbreak size (n) scales as A''^''^. Moreover, the 
maximal and the average duration of an outbreak grow as t* ~ A^^^^ and {t) ~ In A*', respectively. 
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Infection_processes typically involve a threshold 0, 0, 
0, 0, 0, H H Below the epidemic threshold, out- 
breaks quickly die out, while above the threshold, out- 
breaks may take off. We study epidemic outbreaks near 
the threshold. Such outbreaks arise naturally. On the 
one hand, human efforts at disease prevention reduce the 
infection rate thereby crossing the epidemic threshold |2| ■ 
On the other hand, evolution may increase the infection 
rate of diseases hovering just below the threshold, en- 
hancing the likelihood of near-threshold outbreaks [lol| . 
Typically, detection, modeling, and eradication of infec- 
tious diseases are subtle for outbreaks near the epidemic 
threshold. 

The total number of infected individuals is a basic mea- 
sure of the severity of an epidemic outbreak. We study 
outbreaks originating from a single infected individual in 
a large susceptible population. Our main result is that 
near the epidemic threshold, the maximal outbreak size 
grows as a power-law of the population size iV, 
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In contrast, below the epidemic threshold, endemic out- 
breaks involve a small number of infected individuals, 
while above the epidemic threshold, pandemic outbreaks 
involve a fraction of the population ~ N . Therefore, 
outbreaks near the epidemic threshold have a distinct 
intermediate size between a pandemic and an endemic 
outbreak Loosely speaking, epidemics come in three 
sizes: large, medium, and small. 

The scaling law has several important implications 
concerning the statistics of both the size and the duration 
of the outbreaks. It implies that the average size of out- 
breaks (n) and the maximal duration of outbreaks t^, both 
scale as (n) t^, '-^ iV^/^ near the epidemic threshold. 
Furthermore, the average duration of the outbreaks {t) 
scales logarithmically, {t) ~ InA^. These behaviors hold 
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in a sizable range of infection rates, namely in a window 
of the order 0{N~'^/'^) around the epidemic threshold. 

These scaling laws are demonstrated for the classic 
Susceptible-Infected-Recovered (SIR) infection process 
P, 0, In this model, the population consists of s 
susceptible, i infected, and r recovered individuals with 
N = s ~\- i ~\- r. These sub-populations change due to 
two competing processes: infection and recovery. The 
disease is transmitted from an infected individual to a 
susceptible one with rate a/TV, where a is the infection 
rate: 



, . , asi/N . 1 ■ 1 N 

(s,i,r) — y (s - 1,2 + l,r). 



Infected individuals recover with a unit rate: 



(s,z,r) (s,i - l,r+ 1). 
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The infection process starts with a single infected indi- 
vidual, (s,i,r) = (iV — 1,1,0), and it ends with none 
(s, i, r) — {N — n, 0, n). 

The total size of the outbreak n and the duration of 
the outbreak t are the outcomes of a stochastic process. 
We study statistical properties of these random variables, 
particularly, their average and maximal size, as a function 
of the population size (We implicitly consider an average 
over infinitely many realizations of the infection process.) 

In the infinite population limit, the epidemic threshold 
is a = 1. Since infection occurs with probability and 
recovery with probability y^p^ , the average outbreak size 
satisfies (n) = tt — V'^^rr—kn)- Thus, below the thresh- 
old (a < 1), a finite number of individuals is infected, 
(ri) = (1 — a)^^. Above the threshold (a > 1), there is 
a pandemic outbreak with a finite fraction of the pop- 
ulation infected: (n) = rN 0, ^|. At the threshold 
(a = 1), the probability that the outbreak size equals n, 
Gn, is found recursively: Gn = | X]m=i GmGn-m start- 
ing with Gi = 1/2. This recursion reflects the fact 
that the first infection event results in two indepen- 
dent infection processes 0|. The generating function 
is E„>i Gnz^ = l-VT 



from which the size distri- 



2 



A 
c 
V 




N 

FIG. 1: The average outbreak size versus the population 
size for the SIR infection process at the epidemic threshold 
(a = 1). Shown are Monte Carlo simulation results repre- 
senting an average over 10^ independent realizations of the 
infection process (circles). A line of slope 1/3 is also shown 
as a reference. A least-square-fit to (n) ~ A^^ in the range 
10=^ < A < 10^ yields 7 = 0.334 ± 0.001. 



bution is a power-law, 
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for sufficiently large outbreaks n ^ 1. 

For a finite, yet large population, the outbreak size dis- 
tribution Q holds, but only up to the maximal outbreak 
size: 1 <ti n ^ n^. Outbreaks beyond the maximal size 
are practically impossible. Therefore, the average out- 
break size grows according to (n) — J2n=i ~ nl^^. 
Naively assuming that a finite fraction of the population 
may become infected, ~ N, would lead to (n) ~ N^^^. 
While consistent with the generic statistical uncertain- 
ties, this law is in fact erroneous. Instead, the outbreak 
size is much smaller because the epidemic outbreak weak- 
ens as more individuals become infected, and it finally 
dies out when the number of infected individuals becomes 
of the order n* . When there are N ~ n^, susceptible indi- 
viduals, the total infection rate a{N — n^,)i/N shows that 
the infection rate is effectively reduced, acft = I — n^/N. 
Therefore, the epidemic becomes essentially endemic. 
(This is clearly a finite population effect: the suscep- 
tible population "reservoir" is never affected in the in- 
finite population limit.) Equating the outbreak size in 

the endemic phase (n) ^ (1 — acs)'^ ^ N/n^, with that 

1/2 

estimated from the size distribution, (n) ~ n* , gives 
the scaling law ^ governing the maximal outbreak size. 
Hence, in the worse case scenario, only a fraction of the 
order of N~^/^ of the entire population can ever be in- 
fected. 

As a byproduct we obtain the scaling law for the aver- 
age outbreak size 
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FIG. 2: The normalized cumulative distribution 
U„{N) /U„{oo) versus the normalized outbreak size n/N^^^. 
The data corresponds to an average over 10^ independent 
realizations. 



Large scale Monte Carlo simulations confirm this behav- 
ior (Fig. 1). The simulations are a straightforward real- 
ization of the infection process. When there are s sus- 
ceptible individuals, with probability l/{l + as/N) a re- 
covery event occurs, and otherwise, an infection event 
occurs. The simulation results represent an average over 
a remarkably large number of independent realizations. 

Statistical properties of the outbreak size are self- 
similar as they follow a universal, population-size in- 
dependent law. Once the outbreak size distribution 
and the outbreak size are properly normalized by the 
infinite population distribution and the maximal out- 
break size respectively, a universal behavior emerges: 
Gn{N)/Gn{oo) ^ g{Ti/N^/^). This universality, remi- 
niscent of finite-size scaling in critical phenomena i;^], 
was confirmed numerically by studying the cumulative 
distribution UniN) = Em>„G„(7V) (Fig. 2). This pro- 
vides further verification of the scaling law 

The scaling laws characterizing the outbreak size hold 
not only at the threshold but also in a window around 
the threshold. Equating the average outbreak size (O 
with the behavior in the endemic phase, (n) = (1 — a)~^ , 
we find that the threshold window (i.e., the range of in- 
fection rates for which the intermediate behavior holds) 
diminishes with the population size as 
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This parameter range can be sizable for moderate pop- 
ulations — for example, when N — 10'^, the threshold 
window is roughly 0.9 < a < 1.1 and the maximal out- 
break size is smaller than the population size by a factor 
of 10. 

The behavior of (n) near the epidemic thresh- 
old provides another manifestation of the scaling law 
©. Indeed, plotting the average outbreak size ver- 
sus the infection rate normalized according to (O 
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FIG. 3: The near threshold behavior. Shown is the normal- 
ized outbreak size (n) /N^^^ versus the normalized distance 
from the threshold (q — 1)A''^/'^. The data corresponds to an 
average over 10^ independent realizations. 



and ©J respectively, shows a universal behavior: 
(n)/7Vi/3 ^ Q [(1 - a)iVi/3] (Fig. 3). 

The threshold window is larger than the canonical 
N^^/'^ estimate arising either from the standard large- 
population analysis 0, 0] or from the widely-used de- 
terministic SIR ordinary differential equations de- 
scribing the evolution of the average susceptible and in- 
fected populations 0|. Moreover, the related SI (some- 
times also termed SIS) model, where a recovered indi- 
vidual immediately becomes susceptible, is characterized 
by the simpler behavior n* ~ TV and (n) '--^ N^^"^; finite 
size effects are not as pronounced because there is no 
depletion of the susceptible reservoir. 

The scaling laws for the outbreak size have direct im- 
plications concerning the dynamics and in particular, the 
duration of infection processes near the epidemic thresh- 
old. To obtain these scaling laws, we again consider first 
the infinite population limit. At the epidemic threshold, 
a = 1, infection and recovery occur with equal prob- 
abilities and therefore, the average number of infected 
individuals is conserved, I{t) = 1. The probability Pi{t) 
that there are i infected individuals at time t satisfies 



dt 
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together with the initial condition Pi{0) = Sn. The dis- 
tribution is geometric, P,(t) = t'-^{l + HI Ell 
for i > 1, and Po{t) = ^(1-^*)-^ for i = 0. Therefore, the 
survival probability of the outbreak, i.e., the probability 
that the outbreak is still active at time t is simply 



p{t) = (i + ty 



(8) 



since P{t) = 1 — _Po(^)- Restricting attention to active 
outbreaks, the average number of infected individuals 
grows linearly with time (i) = I{t)/P{t) = 1 + t. Con- 
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FIG. 4: The average outbreak duration at the epidemic 
threshold versus the population size. Simulation results, ob- 
tained from an average over 10'' realizations are consistent 
with the theoretical prediction 1101 . A best fit to (t) = /3\nN 
yields l3 = 0.32 ± 0.01. 



sequently, the typical number of recovered individuals 
r ~ dt' (1 -I- 1') grows quadratically with time: r ^ . 

For finite populations, the probability that the out- 
break is still alive at time t decays as P(i, N) ~ up to 
the maximal time scale t <^t^. The survival probability 
is sharply suppressed for times larger than the maximal 
time. The maximal duration of outbreaks is estimated 
by equating the time dependent outbreak size n ^ r ^ 
with the maximal outbreak size ~ iV2/3. Therefore, 



(9) 



The maximal duration of outbreaks greatly exceeds 
both the typical duration that is of the order of one and 
the average duration of an outbreak {t) which exhibits 
an interesting logarithmic growth. To derive the loga- 
rithmic law, we first note that, by definition, the average 
duration of an outbreak is (t) = dtt ■^P{t,N). Us- 
ing the infinite population result ^ and integrating up 
to tit that plays the role of a cutoff, we get 



{t)~-\nN. 



(10) 



Numerical simulations confirm this behavior (Fig. 4). 
The probability distribution for the duration of out- 
breaks also follows a population-size independent law: 
P{t,N)/P{t) (t/N^^^) as shown in Fig. 5. How- 
ever, the convergence to this law is not uniform: it is 
slow for short durations but fast at large durations. 

In summary, we found that outbreaks in the vicinity of 
the epidemic threshold have a distinct size, characterized 
by a distinct power-law dependence of the population 
size. This behavior describes a range of infection pro- 
cesses in the vicinity of the epidemic threshold. The size 
of this threshold window is larger than expected from the 
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FIG. 5: The survival probability at the epidemic threshold. 
Shown is the normalized survival probability P{t,N) / P{t) 
versus the normalized duration time t/N^^^ . The data corre- 
sponds to an average over 10* realizations. 



processes. Typically, one has to compute Pi^n the proba- 
bility that there are i infected individuals and r removed 
individuals from the master equations. Although there 
are N'^ such coupled ordinary differential equations, the 
scaling laws i ^ TV^/'^ and r ^ N^/^ imply that the num- 
ber of relevant equations is much smaller and scales only 
linearly with the population size. 

Several questions arise, e.g., what is the shape of the 
scaling functions Q (n/N^^^) andV (t/N^^'^) characteriz- 
ing the size and duration of outbreaks near the epidemic 
threshold? Analytical determination of these functions is 
very challenging as it requires treatment of the full mas- 
ter equations describing the stochastic infection process 
0, that is, the distribution Pi^r{t,N) is needed [l^ . 

Further related problems include the corresponding 
near-threshold scaling laws for spatial epidemic models, 
where the geometry and the spatial structure of the in- 
fected domain play a role |20, IsJ, |22| , and infection 
processes on networks [23, l24| . 



traditional large system size analysis techniques or from 
the deterministic description. We conclude that statis- 
tical fluctuations and finite population effects are most 
pronounced and may be quite subtle near the epidemic 
threshold. 

The scaling laws have concrete implications regarding 
the computational complexity of near-threshold infection 
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